Single-Cell Profiling Reveals Immune Aberrations in Progressive Idiopathic Pulmonary Fibrosis.

RATIONALE: Changes in peripheral blood cell populations have been observed but not detailed at single-cell resolution in idiopathic pulmonary fibrosis (IPF). OBJECTIVES: To provide an atlas of the changes in the peripheral immune system in stable and progressive IPF. METHODS: Peripheral blood mononuclear cells (PBMCs) from IPF patients and controls were profiled using 10x Chromium 5' single-cell RNA sequencing (scRNA-seq). Flow cytometry was used for validation. Protein concentrations of Regulatory T-cells (Tregs) and Monocytes chemoattractants were measured in plasma and lung homogenates from patients and controls. MEASUREMENTS AND MAIN RESULTS: Thirty-eight PBMC samples from 25 patients with IPF and 13 matched controls yielded 149,564 cells that segregated into 23 subpopulations. Classical monocytes were increased in progressive and stable IPF compared to controls (32.1%, 25.2%, 17.9%, respectively, p<0.05). Total lymphocytes were decreased in IPF vs controls, and in progressive vs stable IPF (52.6% vs 62.6%, p=0.035). Tregs were increased in progressive vs stable IPF (1.8% vs 1.1% of all PBMC, p=0.007), although not different than controls, and may be associated with decreased survival (P=0.009 in Kaplan-Meier analysis; P=0.069 after adjusting for age, sex, and baseline FVC). Flow cytometry analysis confirmed this finding in an independent cohort of IPF patients. Fraction of Tregs out of all T cells was also increased in two cohorts of lung scRNA-seq. CCL22 and CCL18, ligands for CCR4 and CCR8 Treg chemotaxis receptors, were increased in IPF. CONCLUSIONS: The single-cell atlas of the peripheral immune system in IPF, reveals an outcome-predictive increase in classical monocytes and Tregs, as well as evidence for a lung-blood immune recruitment axis involving CCL7 (for classical monocytes) and CCL18/CCL22 (for Tregs).

The regulation and differentiation of regulatory T cells and their dysfunction in autoimmune diseases.

The discovery of FOXP3+ regulatory T (Treg) cells as a distinct cell lineage with a central role in regulating immune responses provided a deeper understanding of self-tolerance. The transcription factor FOXP3 serves a key role in Treg cell lineage determination and maintenance, but is not sufficient to enable the full potential of Treg cell suppression, indicating that other factors orchestrate the fine-tuning of Treg cell function. Moreover, FOXP3-independent mechanisms have recently been shown to contribute to Treg cell dysfunction. FOXP3 mutations in humans cause lethal fulminant systemic autoinflammation (IPEX syndrome). However, it remains unclear to what degree Treg cell dysfunction is contributing to the pathophysiology of common autoimmune diseases. In this Review, we discuss the origins of Treg cells in the periphery and the multilayered mechanisms by which Treg cells are induced, as well as the FOXP3-dependent and FOXP3-independent cellular programmes that maintain the suppressive function of Treg cells in humans and mice. Further, we examine evidence for Treg cell dysfunction in the context of common autoimmune diseases such as multiple sclerosis, inflammatory bowel disease, systemic lupus erythematosus and rheumatoid arthritis.

Systems Analysis of Immune Changes after B-cell Depletion in Autoimmune Multiple Sclerosis.

Multiple sclerosis (MS) is a complex genetically mediated autoimmune disease of the central nervous system where anti-CD20-mediated B cell depletion is remarkably effective in the treatment of early disease. While previous studies investigated the effect of B cell depletion on select immune cell subsets using flow cytometry-based methods, the therapeutic impact on patient immune landscape is unknown. In this study, we explored how a therapy-driven " in vivo perturbation " modulates the diverse immune landscape by measuring transcriptomic granularity with single-cell RNA sequencing (scRNAseq). We demonstrate that B cell depletion leads to cell type-specific changes in the abundance and function of CSF macrophages and peripheral blood monocytes. Specifically, a CSF-specific macrophage population with an anti-inflammatory transcriptomic signature and peripheral CD16 + monocytes increased in frequency post-B cell depletion. In addition, we observed increases in TNFα messenger RNA and protein in monocytes post-B cell depletion, consistent with the finding that anti-TNFα treatment exacerbates autoimmune activity in MS. In parallel, B cell depletion also induced changes in peripheral CD4 + T cell populations, including increases in the frequency of TIGIT + regulatory T cells and marked decreases in the frequency of myelin peptide loaded-tetramer binding CD4 + T cells. Collectively, this study provides an exhaustive transcriptomic map of immunological changes, revealing different mechanisms of action contributing to the high efficacy in B cell depletion treatment of MS.

Peripheral Blood Single-Cell Sequencing Uncovers Common and Specific Immune Aberrations in Fibrotic Lung Diseases.

Rationale and Objectives: The extent and commonality of peripheral blood immune aberrations in fibrotic interstitial lung diseases are not well characterized. In this study, we aimed to identify common and distinct immune aberrations in patients with idiopathic pulmonary fibrosis (IPF) and fibrotic hypersensitivity pneumonitis (FHP) using cutting-edge single-cell profiling technologies. Methods: Single-cell RNA sequencing was performed on patients and healthy controls' peripheral blood and bronchoalveolar lavage samples using 10X Genomics 5' gene expression and V(D)J profiling. Cell type composition, transcriptional profiles, cellular trajectories and signaling, and T and B cell receptor repertoires were studied. The standard Seurat R pipeline was followed for cell type composition and differential gene expression analyses. Transcription factor activity was imputed using the DoRothEA-VIPER algorithm. Pseudotime analyses were conducted using Monocle3, while RNA velocity analyses were performed with Velocyto, scVelo, and CellRank. Cell-cell connectomics were assessed using the Connectome R package. V(D)J analyses were conducted using CellRanger and Immcantation frameworks. Across all analyses, disease group differences were assessed using the Wilcoxon rank-sum test. Measurements and Main Results: 327,990 cells from 83 samples were profiled. Overall, changes in monocytes were common to IPF and FHP, whereas lymphocytes exhibited disease-specific aberrations. Both diseases displayed enrichment of CCL3 hi /CCL4 hi CD14+ monocytes (p<2.2e-16) and S100A hi CD14+ monocytes (p<2.2e-16) versus controls. Trajectory and RNA velocity analysis suggested that pro-fibrotic macrophages observed in BAL originated from peripheral blood monocytes. Lymphocytes exhibited disease-specific aberrations, with CD8+ GZMK hi T cells and activated B cells primarily enriched in FHP patients. V(D)J analyses revealed unique T and B cell receptor complementarity-determining region 3 (CDR3) amino acid compositions (p<0.05) in FHP and significant IgA enrichment in IPF (p<5.2e-7). Conclusions: We identified common and disease-specific immune mechanisms in IPF and FHP; S100A hi monocytes and SPP1 hi macrophages are common to IPF and FHP, whereas GMZK hi T lymphocytes and T and B cell receptor repertoires were unique in FHP. Our findings open novel strategies for the diagnosis and treatment of IPF and FHP.

Latent EBV impairs immune cell signaling and enhances the efficacy of anti-CD3 mAb in Type 1 Diabetes.

Teplizumab has been approved for the delay of the onset of type 1 diabetes and may modulate new onset disease. We found that patients who were EBV positive at baseline had a more robust response to drug in two clinical trials and therefore postulated that latent virus has general effects in modifying immune responses. We compared the phenotypes, transcriptomes, and development of peripheral blood cells before and after teplizumab treatment. Higher number of Tregs and partially exhausted CD8 + T cells were found in EBV seropositive individuals at the baseline in the TN10 trial and AbATE trial. Single cell transcriptomics and functional assays identified downregulation of the T cell receptor and other signaling pathways before treatment. Impairments in function of adaptive immune cells were enhanced by teplizumab treatment in EBV seropositive individuals. Our data indicate that EBV can impair signaling pathways generally in immune cells, that broadly redirect cell differentiation.

Cytokinopathy with aberrant cytotoxic lymphocytes and profibrotic myeloid response in SARS-CoV-2 mRNA vaccine-associated myocarditis.

Rare immune-mediated cardiac tissue inflammation can occur after vaccination, including after SARS-CoV-2 mRNA vaccines. However, the underlying immune cellular and molecular mechanisms driving this pathology remain poorly understood. Here, we investigated a cohort of patients who developed myocarditis and/or pericarditis with elevated troponin, B-type natriuretic peptide, and C-reactive protein levels as well as cardiac imaging abnormalities shortly after SARS-CoV-2 mRNA vaccination. Contrary to early hypotheses, patients did not demonstrate features of hypersensitivity myocarditis, nor did they have exaggerated SARS-CoV-2-specific or neutralizing antibody responses consistent with a hyperimmune humoral mechanism. We additionally found no evidence of cardiac-targeted autoantibodies. Instead, unbiased systematic immune serum profiling revealed elevations in circulating interleukins (IL-1ß, IL-1RA, and IL-15), chemokines (CCL4, CXCL1, and CXCL10), and matrix metalloproteases (MMP1, MMP8, MMP9, and TIMP1). Subsequent deep immune profiling using single-cell RNA and repertoire sequencing of peripheral blood mononuclear cells during acute disease revealed expansion of activated CXCR3+ cytotoxic T cells and NK cells, both phenotypically resembling cytokine-driven killer cells. In addition, patients displayed signatures of inflammatory and profibrotic CCR2+ CD163+ monocytes, coupled with elevated serum-soluble CD163, that may be linked to the late gadolinium enhancement on cardiac MRI, which can persist for months after vaccination. Together, our results demonstrate up-regulation in inflammatory cytokines and corresponding lymphocytes with tissue-damaging capabilities, suggesting a cytokine-dependent pathology, which may further be accompanied by myeloid cell-associated cardiac fibrosis. These findings likely rule out some previously proposed mechanisms of mRNA vaccine--associated myopericarditis and point to new ones with relevance to vaccine development and clinical care.

Single-cell profiling reveals immune aberrations in progressive idiopathic pulmonary fibrosis.

Rationale: Changes in peripheral blood cell populations have been observed but not detailed at single-cell resolution in idiopathic pulmonary fibrosis (IPF). Objectives: To provide an atlas of the changes in the peripheral immune system in stable and progressive IPF. Methods: Peripheral blood mononuclear cells (PBMCs) from IPF patients and controls were profiled using 10x Chromium 5' single-cell RNA sequencing (scRNA-seq). Flow cytometry was used for validation. Protein concentrations of Regulatory T-cells (Tregs) and Monocytes chemoattractants were measured in plasma and lung homogenates from patients and controls. Measurements and Main Results: Thirty-eight PBMC samples from 25 patients with IPF and 13 matched controls yielded 149,564 cells that segregated into 23 subpopulations, corresponding to all expected peripheral blood cell populations. Classical monocytes were increased in progressive and stable IPF compared to controls (32.1%, 25.2%, 17.9%, respectively, p<0.05). Total lymphocytes were decreased in IPF vs controls, and in progressive vs stable IPF (52.6% vs 62.6%, p=0.035). Tregs were increased in progressive IPF (1.8% vs 1.1%, p=0.007), and were associated with decreased survival (P=0.009 in Kaplan-Meier analysis). Flow cytometry analysis confirmed this finding in an independent cohort of IPF patients. Tregs were also increased in two cohorts of lung scRNA-seq. CCL22 and CCL18, ligands for CCR4 and CCR8 Treg chemotaxis receptors, were increased in IPF. Conclusions: The single-cell atlas of the peripheral immune system in IPF, reveals an outcome-predictive increase in classical monocytes and Tregs, as well as evidence for a lung-blood immune recruitment axis involving CCL7 (for classical monocytes) and CCL18/CCL22 (for Tregs).

Regulatory T cells in peripheral tissue tolerance and diseases.

Maintenance of peripheral tolerance by CD4+Foxp3+ regulatory T cells (Tregs) is essential for regulating autoreactive T cells. The loss of function of Foxp3 leads to autoimmune disease in both animals and humans. An example is the rare, X-linked recessive disorder known as IPEX (Immune Dysregulation, Polyendocrinopathy, Enteropathy X-linked) syndrome. In more common human autoimmune diseases, defects in Treg function are accompanied with aberrant effector cytokines such as IFNγ. It has recently become appreciated that Tregs plays an important role in not only maintaining immune homeostasis but also in establishing the tissue microenvironment and homeostasis of non-lymphoid tissues. Tissue resident Tregs show profiles that are unique to their local environments which are composed of both immune and non-immune cells. Core tissue-residence gene signatures are shared across different tissue Tregs and are crucial to homeostatic regulation and maintaining the tissue Treg pool in a steady state. Through interaction with immunocytes and non-immunocytes, tissue Tregs exert a suppressive function via conventional ways involving contact dependent and independent processes. In addition, tissue resident Tregs communicate with other tissue resident cells which allows Tregs to adopt to their local microenvironment. These bidirectional interactions are dependent on the specific tissue environment. Here, we summarize the recent advancements of tissue Treg studies in both human and mice, and discuss the molecular mechanisms that maintain tissue homeostasis and prevent pathogenesis.

Hyperosmotic stress response regulates interstitial homeostasis and pathogenic inflammation.

Hyperosmotic stress triggers an evolutionally preserved, fundamental cellular response. A growing body of evidence has highlighted the role of extra-renal, interstitial hyperosmolality in maintaining local tissue immune homeostasis and potentially driving tissue inflammation in human diseases. The hyperosmotic stress response initiates cellular shrinkage, oxidative stress, metabolic remodeling and cell cycle arrest, all of which are adjusted by a counteractive adaptative response that includes osmolyte synthesis, upregulation of ion transporters and induction of heat shock proteins. Recent studies have revealed that high osmolality can impact immune cell differentiation and activation pathways in a cell type specific manner. The fine-tuning of the immune response depends on the tissue microenvironment. Accordingly, novel therapeutic approaches that target hyperosmolality-mediated inflammation may be identified by furthering our understanding of hyperosmotic response in the context of disease. In this review, we discuss the cellular and molecular mechanisms by which hyperosmotic stress response regulates interstitial homeostasis and pathogenic inflammation.

PD-1highCXCR5-CD4+ peripheral helper T cells promote CXCR3+ plasmablasts in human acute viral infection.

T cell-B cell interaction is the key immune response to protect the host from severe viral infection. However, how T cells support B cells to exert protective humoral immunity in humans is not well understood. Here, we use COVID-19 as a model of acute viral infections and analyze CD4+ T cell subsets associated with plasmablast expansion and clinical outcome. Peripheral helper T cells (Tph cells; denoted as PD-1highCXCR5-CD4+ T cells) are significantly increased, as are plasmablasts. Tph cells exhibit "B cell help" signatures and induce plasmablast differentiation in vitro. Interestingly, expanded plasmablasts show increased CXCR3 expression, which is positively correlated with higher frequency of activated Tph cells and better clinical outcome. Mechanistically, Tph cells help B cell differentiation and produce more interferon γ (IFNγ), which induces CXCR3 expression on plasmablasts. These results elucidate a role for Tph cells in regulating protective B cell response during acute viral infection.

A multiple sclerosis-protective coding variant reveals an essential role for HDAC7 in regulatory T cells.

Genome-wide association studies identifying hundreds of susceptibility loci for autoimmune diseases indicate that genes active in immune cells predominantly mediate risk. However, identification and functional characterization of causal variants remain challenging. Here, we focused on the immunomodulatory role of a protective variant of histone deacetylase 7 (HDAC7). This variant (rs148755202, HDAC7.p.R166H) was identified in a study of low-frequency coding variation in multiple sclerosis (MS). Through transcriptomic analyses, we demonstrate that wild-type HDAC7 regulates genes essential for the function of Foxp3+ regulatory T cells (Tregs), an immunosuppressive subset of CD4 T cells that is generally dysfunctional in patients with MS. Moreover, Treg-specific conditional hemizygous deletion of HDAC7 increased the severity of experimental autoimmune encephalitis (EAE), a mouse model of neuroinflammation. In contrast, Tregs transduced with the protective HDAC7 R166H variant exhibited higher suppressive capacity in an in vitro functional assay, mirroring phenotypes previously observed in patient samples. In vivo modeling of the human HDAC7 R166H variant by generation of a knock-in mouse model bearing an orthologous R150H substitution demonstrated decreased EAE severity linked to transcriptomic alterations of brain-infiltrating Tregs, as assessed by single-cell RNA sequencing. Our data suggest that dysregulation of epigenetic modifiers, a distinct molecular class associated with disease risk, may influence disease onset. Last, our approach provides a template for the translation of genetic susceptibility loci to detailed functional characterization, using in vitro and in vivo modeling.

Impaired TIGIT expression on B cells drives circulating follicular helper T cell expansion in multiple sclerosis.

B cell depletion in patients with relapsing-remitting multiple sclerosis (RRMS) markedly prevents new MRI-detected lesions and disease activity, suggesting the hypothesis that altered B cell function leads to the activation of T cells driving disease pathogenesis. Here, we performed comprehensive analyses of CD40 ligand- (CD40L-) and IL-21-stimulated memory B cells from patients with MS and healthy age-matched controls, modeling the help of follicular helper T cells (Tfh cells), and found a differential gene expression signature in multiple B cell pathways. Most striking was the impaired TIGIT expression on MS-derived B cells mediated by dysregulation of the transcription factor TCF4. Activated circulating Tfh cells (cTfh cells) expressed CD155, the ligand of TIGIT, and TIGIT on B cells revealed their capacity to suppress the proliferation of IL-17-producing cTfh cells via the TIGIT/CD155 axis. Finally, CCR6+ cTfh cells were significantly increased in patients with MS, and their frequency was inversely correlated with that of TIGIT+ B cells. Together, these data suggest that the dysregulation of negative feedback loops between TIGIT+ memory B cells and cTfh cells in MS drives the activated immune system in this disease.

Identity thieves: T cells steal CD20 from B cells but mark themselves for certain death.

T cells can acquire CD20 from B cells via trogocytosis, then contribute to autoimmune neurological disease while becoming targets for therapeutically effective B cell depletion.

Population genetics meets single-cell sequencing.

Single-cell technology can be used to understand the genetic basis of human diseases.

Type I interferon transcriptional network regulates expression of coinhibitory receptors in human T cells.

Although inhibition of T cell coinhibitory receptors has revolutionized cancer therapy, the mechanisms governing their expression on human T cells have not been elucidated. In the present study, we show that type 1 interferon (IFN-I) regulates coinhibitory receptor expression on human T cells, inducing PD-1/TIM-3/LAG-3 while inhibiting TIGIT expression. High-temporal-resolution mRNA profiling of IFN-I responses established the dynamic regulatory networks uncovering three temporal transcriptional waves. Perturbation of key transcription factors (TFs) and TF footprint analysis revealed two regulator modules with different temporal kinetics that control expression of coinhibitory receptors and IFN-I response genes, with SP140 highlighted as one of the key regulators that differentiates LAG-3 and TIGIT expression. Finally, we found that the dynamic IFN-I response in vitro closely mirrored T cell features in acute SARS-CoV-2 infection. The identification of unique TFs controlling coinhibitory receptor expression under IFN-I response may provide targets for enhancement of immunotherapy in cancer, infectious diseases and autoimmunity.

Single-cell multi-omics reveals dyssynchrony of the innate and adaptive immune system in progressive COVID-19.

Dysregulated immune responses against the SARS-CoV-2 virus are instrumental in severe COVID-19. However, the immune signatures associated with immunopathology are poorly understood. Here we use multi-omics single-cell analysis to probe the dynamic immune responses in hospitalized patients with stable or progressive course of COVID-19, explore V(D)J repertoires, and assess the cellular effects of tocilizumab. Coordinated profiling of gene expression and cell lineage protein markers shows that S100Ahi/HLA-DRlo classical monocytes and activated LAG-3hi T cells are hallmarks of progressive disease and highlights the abnormal MHC-II/LAG-3 interaction on myeloid and T cells, respectively. We also find skewed T cell receptor repertories in expanded effector CD8+ clones, unmutated IGHG+ B cell clones, and mutated B cell clones with stable somatic hypermutation frequency over time. In conclusion, our in-depth immune profiling reveals dyssynchrony of the innate and adaptive immune interaction in progressive COVID-19.

NEBULA is a fast negative binomial mixed model for differential or co-expression analysis of large-scale multi-subject single-cell data.

The increasing availability of single-cell data revolutionizes the understanding of biological mechanisms at cellular resolution. For differential expression analysis in multi-subject single-cell data, negative binomial mixed models account for both subject-level and cell-level overdispersions, but are computationally demanding. Here, we propose an efficient NEgative Binomial mixed model Using a Large-sample Approximation (NEBULA). The speed gain is achieved by analytically solving high-dimensional integrals instead of using the Laplace approximation. We demonstrate that NEBULA is orders of magnitude faster than existing tools and controls false-positive errors in marker gene identification and co-expression analysis. Using NEBULA in Alzheimer's disease cohort data sets, we found that the cell-level expression of APOE correlated with that of other genetic risk factors (including CLU, CST3, TREM2, C1q, and ITM2B) in a cell-type-specific pattern and an isoform-dependent manner in microglia. NEBULA opens up a new avenue for the broad application of mixed models to large-scale multi-subject single-cell data.

Cutting Edge: Distinct B Cell Repertoires Characterize Patients with Mild and Severe COVID-19.

Protective immunity against COVID-19 likely depends on the production of SARS-CoV-2-specific plasma cells and memory B cells postinfection or postvaccination. Previous work has found that germinal center reactions are disrupted in severe COVID-19. This may adversely affect long-term immunity against reinfection. Consistent with an extrafollicular B cell response, patients with severe COVID-19 have elevated frequencies of clonally expanded, class-switched, unmutated plasmablasts. However, it is unclear whether B cell populations in individuals with mild COVID-19 are similarly skewed. In this study, we use single-cell RNA sequencing of B cells to show that in contrast to patients with severe COVID-19, subjects with mildly symptomatic COVID-19 have B cell repertoires enriched for clonally diverse, somatically hypermutated memory B cells ∼30 d after the onset of symptoms. This provides evidence that B cell responses are less disrupted in mild COVID-19 and result in the production of memory B cells.

Immune dysregulation and autoreactivity correlate with disease severity in SARS-CoV-2-associated multisystem inflammatory syndrome in children.

Multisystem inflammatory syndrome in children (MIS-C) is a life-threatening post-infectious complication occurring unpredictably weeks after mild or asymptomatic SARS-CoV-2 infection. We profiled MIS-C, adult COVID-19, and healthy pediatric and adult individuals using single-cell RNA sequencing, flow cytometry, antigen receptor repertoire analysis, and unbiased serum proteomics, which collectively identified a signature in MIS-C patients that correlated with disease severity. Despite having no evidence of active infection, MIS-C patients had elevated S100A-family alarmins and decreased antigen presentation signatures, indicative of myeloid dysfunction. MIS-C patients showed elevated expression of cytotoxicity genes in NK and CD8+ T cells and expansion of specific IgG-expressing plasmablasts. Clinically severe MIS-C patients displayed skewed memory T cell TCR repertoires and autoimmunity characterized by endothelium-reactive IgG. The alarmin, cytotoxicity, TCR repertoire, and plasmablast signatures we defined have potential for application in the clinic to better diagnose and potentially predict disease severity early in the course of MIS-C.

Post-infectious inflammatory disease in MIS-C features elevated cytotoxicity signatures and autoreactivity that correlates with severity.

Multisystem inflammatory syndrome in children (MIS-C) is a life-threatening post-infectious complication occurring unpredictably weeks after mild or asymptomatic SARS-CoV2 infection in otherwise healthy children. Here, we define immune abnormalities in MIS-C compared to adult COVID-19 and pediatric/adult healthy controls using single-cell RNA sequencing, antigen receptor repertoire analysis, unbiased serum proteomics, and in vitro assays. Despite no evidence of active infection, we uncover elevated S100A-family alarmins in myeloid cells and marked enrichment of serum proteins that map to myeloid cells and pathways including cytokines, complement/coagulation, and fluid shear stress in MIS-C patients. Moreover, NK and CD8 T cell cytotoxicity genes are elevated, and plasmablasts harboring IgG1 and IgG3 are expanded. Consistently, we detect elevated binding of serum IgG from severe MIS-C patients to activated human cardiac microvascular endothelial cells in culture. Thus, we define immunopathology features of MIS-C with implications for predicting and managing this SARS-CoV2-induced critical illness in children.